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Microrefrigerators based on normal metal-insulator-superconductor (NIS) junctions represent a 
very attractive alternative to cool the microbolometers and calorimeters for astrophysical observa- 
tions in space-borne experiments. The performance in such measurements requires a good knowledge 
of the noise sources in the detectors. In this paper we present detailed calculations of the thermal 
fluctuations and of the noise equivalent power due to the heat transfer through the NIS junctions 
or due to the thermal contact between different subsystems of the detector. The influence of the 
background radiation will also be evaluated. Analytical approximations, valid at low temperatures, 
are given. 

I. INTRODUCTION 

The most sensitive detectors, suitable for observation of cosmic sources in far-infrared and X-ray bands, are the 
cryogenic bolometers and microcalorimeters. These detectors operate typically at a temperature of about 100 mK. 
An attractive method for the last stage of detector cooling (from 300 mK to 100 mK) in space-borne experiments is 
based on normal metal-insulator-superconductor (NIS) refrigerating junctions The sensitivity required by such 

experiments motivates the study of the noise sources in the detectors. 

A thermal detector system, such as a bolometer or a calorimeter, consists of a thermal sensing element (TSE) which 
is connected to a heat sink. The TSE typically consists of an absorber and a thermometer. The thermometer can 
be a transition-edge sensor or a (SI)NIS tunnel junction thermometer [^Q. In the rest of this paper we shall 

study the temperature fluctuations and the noise equivalent power {NEP) in the electron system (ES) - which is 
the actual sensing device - of the TSE. In the presence of the optical load due to the background radiation, in the 
concrete calculations, the aim will be to keep the temperature of the ES, Te, around 0.1 K. The temperature of the 
heat bath, T2, could be significantly higher then To, say T2 ~ 0.3 K, in which case the refrigerating junctions are 
connected directly to the TSE (the direct cooling method), or T2 could be about the same as Tg, say T2 ~ To ~ 0.1 K, 
when a refrigerating system is used to cool the heat bath (the indirect cooling method) and, eventually, another 
set of junctions connected to the TSE is used to compensate for the optical load (this will be called the combined 
cooling method). 

Let us start by supposing that a set of NIS cooling junctions is connected to the TSE and calculate the tunneling 
currents, the power extracted from the ES of the TSE, and the shot noise in the tunneling current. Following Ref. 
pof we define four tunneling currents between the ES of the normal TSE and the superconductor (we assume that 
the tunnel barrier is sufficiently thick, so we neglect the Andreev reflection of electrons between the normal metal and 
the superconductor |jll|): 

ji(e)=g(e)/o(e- el/,To)[l - Me,T,)]/e^RT, 
j2{e)=9ie)Me + eV,T,)[l - Me,T,)]/e^RT, 

J3(e)=ff(e)[l - /o(e- eV-,To)]/,(e,T,)/e2i?T, ^ ' 

Me)=g{e)[l - Me + eV,%)]U{e,T,)/e^R^, 



where /e,s(e, Te,s) represent the populations of the electron (in the normal metal) and quasiparticle (in the superconduc- 
tor) energy levels (T, is an effective temperature in the superconductor, used to describe the population of the quasipar- 
ticle energy levels). V is the voltage across the junction, e is the elementary charge, and 5(e) = 9{\e\ — A)|e|/Ve^ — 
is the normalized density of states in the superconductor, where is the Heaviside step function. The energy e is mea- 
sured from the Fermi energy in the superconductor and in Eqs. (|l]) it is always taken in absolute value. Although /e_s 
may not be Fermi distributions in our case of nonequilibrium |12| ] we make the assumption that 1— /o(— e. To) = /e(e, Te) 
(which is an identity for a Fermi distribution) to transform the expressions that involved negative e. In what follows 
we shall concentrate on the junction where the flux of electrons is oriented from the normal metal into the supercon- 
ductor, where eV is positive. Using the definitions given in Eqs. (El), the particle and excitation fluxes, iVj and A^'j, 
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respectively, can be written as 

= -/c - / (ii - 32 - J3 + ]i) de (2) 



e 

1 



;(e)[/c(e " eV, T,) ~ Me + eV, T,)] de, 



^'j = / {n+]2-n~k)de (3) 



A 

1 



.9(e) [/e(e - eV, %) + /^(e + eV, T,) - 2/,(e, %)] de. 



If we denote by ep the Fermi energy of the ES in the TSE, then the total power extracted through a NIS junction 
can be written as 

[(e - eV + e^)[ji - J3) + {ep - e - eV){j4 - J2)] de (4) 

A 

= Qj+epN], 

where 

/■OO 

Qj = / [{e-eV)iji-j3)~{e + eV}{j4-j2)]de (5) 

J A 

00 

+ j2 - js - ji) - eV{ji - j2 - js + 34)] de. 



A 

Since the coolers consist usually of pairs of NIS junctions biased in opposite directions [Q, the total power should be 
calculated as a sum of terms like the ones given in Eq. (^. In a stationary case the total number of electrons in the 
TSE should be constant and the terms of the type epNj cancel each other and we are left with the sum over Qj's. 

The spectral density of the shot noise in the power extracted through the NIS junction [given by Eq. (^)] can be 
written as: 



(S'Q'j)^ = 2 / [{e-eV + eFy{n+h) + ieF-e-eVy{ji+j2)]de 
deg{e)l{e-eV)' 



X [/(e - eV, r,)(l - 2/(e, T,)) + f{e, T,)] 

+ (6 + eV flfie + eV, Te)(l - 2/(e, T,)) + /(e, T,)] [ (6) 
2 



deg{e) 



^JA 

+ eV, Te)(l - 2/(e, T,)) + 2/(e, T,) 

4 



/(e-e\/,Te)(l-2/(6,T,)) 

(7) 



ep deig{e)ie-eV) 



c^Rt' 

X [fie - eV, re)(l - 2/(e, T,)) + /(e, T,)] 

-(e + eV)[f{e + eV, T,){1 - 2/(e, T,)) + /(e, T,)]| (8) 

(see Chap. 1, for a general introduction). In the construction of microbolometers the TSE is sufficiently large 
and the Coulomb blockade phenomenon cannot influence the tunneling of individual electrons. In such a case, in 
the calculation of the spectral density of the shot noise in the total power flux, the contributions of the terms (|^), 
(0), and (H), corresponding to different junctions are simply added to each other. As a result, one would obtain a 
fluctuation in power which is many orders of magnitude larger than what was previously calculated p^ , p^ and which 
corresponds just to the term labeled as Eq. Nevertheless, we shall show next that the extra terms[(0) and (||)] 
have no effect on the observable quantities. The shot noise in the output power has no direct connection with the 
NEP of the detector, but the fluctuations in the power would produce fluctuations in the electronic temperature. 
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which depreciate the detection properties of the microbolometer. The NEP is then defined as the minimum power 
of the input signal that would produce a change in temperature equal to the square root of the mean square fluctuation 
of temperature. 



II. NOISE IN THE THERMAL SENSING ELEMENT 



A. Temperature fluctuations 

The detection of radiation using microbolometers is based on the measurement of quantities which depend on the 
electronic temperature in the detector and, eventually, on the chemical potential of the electrons (as is the case of 
NIS junctions used as thermometers). Therefore, the noise in the total electronic energy is not directly connected to 
the figure of merit of the detector, which is the NEP. 

At very low temperatures, the chemical potential of a Fermi system is related to the Fermi energy by the relation 
/i w ep [l — (7r^/12)(/cBT/eF)^] , where we used some obvious notations. On the other hand, the Fermi energy is 

determined by the density of particles, n = N/V, in the following way: n = {2m^/h^f/^el^'^/{3Tr^). Other quantities 
of interest are the energy per particle u and the specific heat at constant volume, cv = [du/ dT) ^constant, which, 
in the low temperature hmit have the expressions u = (3/5)eF[l + (57r^/12)(A:Br/eF)^] and cy = (7r^/2)(fc|T/eF). 

Let us now calculate the change in the temperature of the electron gas when SNf: electrons of energy e + cp leave 
the normal metal. This loss of electrons changes u by an amount: 

^ U~SN,{e + eF) U 6N, , ^SN, 
'^"^ N-SN, "iv""— '^^ + ^^) — ' 

which can be expressed as 

The change in the Fermi energy can be written as Sep = —{2/3)eF6N/N. Using the two equations we have for Su, 
[Eqs. and (^0|)], and neglecting the terms of the order a"^ = (k-QT/ep)'^ <C 1, we end up with an expression for 
the temperature variation: 

Cy6T, = -J-^. (11) 

Here the subscript e in the notation dT^ refers to the fact that the variation in temperature is due to the loss of 
electrons from the level with energy e + ep. In the linear approximation, the temperature change is an additive 
quantity. Taking the time derivative on both sides of Eq. (|l^) , we can write: 

CySt = -e^. (12) 



Using the expressions for the chemical potential, the internal energy and Eq. (^Tj) , it can be shown that the variation 
of /i with temperature introduces corrections of order a"^ to the first order terms, in all the following equations. 
Therefore in the rest of this paper we shall consider fj. ^ ep. 



B. Noise Equivalent Power 

Having now calculated the variations in temperature due to the variations in the population of the energy levels, 
we can determine the total temperature fluctuation in the ES of the TSE. The system is represented schematically 
in Fig. 1^. The electrons interact with the lattice, at temperature Ti, the detector lattice exchanges heat with a heat 
bath, at temperature T2 , and the thermal resistance between the lattice and the heat bath is the Kapitza resistance. 
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FIG. 1. A schematic drawing of the thermal sensing element. The two rectangular boxes represent two distinct subsystems 
of this element: the electron system and the lattice. The electron system, the lattice, the heat bath, and the superconductor 
are at the temperatures Te, Ti, T2, and Tg, respectively. In the case of direct cooling the superconductor is connected with 
the electron gas through NIS junctions, while in the case of indirect cooling it is used to cool down the thermal bath. The 
heat capacities of the electron gas and of the lattice are Cvo and Cvi , respectively. The power fluxes transmitted between the 
systems are represented in the figure by long arrows and are denoted by Qj, Qcp, Qk- Qb is the bias power, while Qoe and 
Qop are the optical input powers into the electron system and into the lattice, respectively. 



We have to explain here the notation. Q will be used for the power fluxes (with the directions given by the long 
arrows in Fig. |l|) and Ng for the time variation of the occupation number of the electronic level of energy e (taken 
with changed sign sign, as we did in the previous section). The subscripts are used to specify the processes that 
caused the power transfer or the time variation of the particle population. In this way, J refers to the particle transfer 
through the NIS junction, ep to the electron-phonon interaction in the TSE, K to the heat transfer through the 
contact between the lattice of the TSE and the heat bath, while oe and op refer to the optical input power, due to the 
external radiation, into the ES and lattice of the TSE, respectively. By Qb we denote the bias power, as shown in the 
figure, but this contribution will be disregarded until the end of the paper, where it will be discussed in connection 
with the uncertainty introduced in experiment by the measurement process. Moreover, we shall use the Greek letter 
i5 in front of the notation to specify the fluctuations of a quantity and the subscript shot for the fluctuations of this 
quantity, due solely to the shot noise (discrete amounts transported randomly at a constant average rate) . 

If we add up all the noise terms into the power balance equation, disregarding the fluctuations of T2 and Tg, we 
arrive at the following set of equations: 




6N, 



9(Ar,j+jV.,p+jV,„J 



^cp ,sliot 



9T„ 



dN, 



oc,Bhot J 



?Gp,shot + <5Qk 



shot 



^Qop,shot 



dTi 



-STi 



(13) 



9(Qep~QK-Qop) 

dTi 



ST, 



■ST, 



where, transforming the summations over the energy levels into integrals, we use the notation ST, = J deaoST^^ in 
which ctq is the density of the energy levels of the electrons in the normal metal, at Fermi energy (we make here the 
usual assumption that the energy range of interest is much smaller than the Fermi energy, so we can consider the 
density of states as being almost constant) . Calculating the Fourier transformations of the two equations above and 
replacing the expression for STi(uj) obtained from the second equation into the flrst equation, we get 



iLL!CycST^^{(^) + STc{uj) 



an 



dTi OTe 



2ajCvi + 



£(Qcp__£k-Qop) 



4 



, SNt^p -l5Qap,Bhot ('^) + SQK,shat ) +<5Qop .shot (l^) 

- OTi 



We introduce the notation: 



iojOviH grj 



We can now integrate both sides of Eq. (|l5|) and write 

d ,^ 



iujCyc + -^{Qj - Qcp - Qc 



iuCyi H qtjt- 

Integrating the r. h. s. of Eq. ( p^ we obtain another expression for At{lli): 

At{uj) = -(5(3j,shot(w) + ^(3op,shot(w) + 5Qoc,shot(w) 

dep dQj '5iVj,shot(^) 
dN dE iiu 

(^) 

^ zwovi H 37^;^ 

Using Eqs. (|l|) and we obtain the foUowing expressions for the mean square value of At(w): 



1 dQcp dQcp 

3(Qcp-Qk-Qop) dTc dTi 



iujCye + ~^^{Qi ~ Qcp — Qoc) 



dQcp dQep 



9Ti 



('^■^Qcp^shot) 



((5 (3K,shot)c^ + (<5 Qop,shot) 



dQcp 
dTi 



+ {6^Q oc,shot/L 



d{Qcp ~Qk~Q 



op) 



^Qj.shot + 



dep dQ] SNj, 



,shot 



dN dE iuj 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



The last term of Eq. (^9|), say T(a;), can be calculated. The correlation function between (5(3j,shot(t^) and 
SNj,shot{ijj)/iijJ gives no contribution due to the 7r/2 phase difference, and the final result is: 



,shot /a; 



(20) 



From the expressions ( p^ ) and (|T^) , using the result (^0|) , one can calculate the spectral density of the temperature 



noise. 
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To determine the amplitude of the signal power that has to be introduced into the detector, to produce changes in 
the temperature of the ES equal to the fluctuations calculated above, we write a set of equations similar to (|l3[): 

CycSTc — Qs — QjO + QcpO + QocO (21) 

d{Qj — Qcp — Qoc) ™ , dQcp ^ 

m ^^'^ + -m^^' ' 

r AT ~ ri h m ^(Qk + Qopo - Qep) 

OviOJl — VKO + LtfopO ~ ^opO H 

Oil 

-^-^Te. (22) 

Here Qs is the power of the input signal into the TSE. The subscript is added to the usual subscripts to denote 
the equilibrium values of the powers defined in Fig. 0. Therefore, — Qjo + Qopo + QocO Qko + Qopo — Qcvo — 0- 
Calculating the Fourier transformations of Eqs. ( pl| ) and (|2|) and replacing STi{u!) from Eq. ( p2[ ) into Eq. (pTl), we 
obtain the Noise Equivalent Power, as defined in the end of Section |: 

NEP^ = \Q,{uj)\^ (23) 

Therefore, NEP can be calculated directly from Eq. (19), with the use of Eq. ([20|). In this way we observe that 
when the detector is cooled indirectly (Qj = and T2 sa Ti w To), T = and NEP is lower. We shall evaluate the 
total NEP in the next sections. 



C. Calculation of the shot noise terms 

a. Electron-phonon shot noise. The heat power transferred between the ES and the lattice of the TSE, Qcp, due 
to the electron-phonon interaction, can be calculated using the formula 

Qep = Sepr!(T5-T5), (24) 

where f2 is the volume of the TSE and Ecp is the electron-phonon coupling constant p6[ |. If, for example, the TSE is 
made of copper, then Sop ~ 4 nW/K^ /im^ 0|. Using the same model as in Ref. the shot noise of Qop has been 
evaluated in Ref. |]l^ . An approximative expression for this power shot noise, namely 

{S^Qcp,shot)u ~ 5Sopr!(T,^ + ), (25) 

has been used in Ref. fill (although a factor of two was introduced by mistake there) . The formula above is identical 
to the exact expression for To = Ti and deviates from it with less than 2% for any To < Ti. Therefore we shall use 
Eq. (p5| ) in the rest of this paper, since its accuracy is good enough for our purposes. 

b. Kapitza shot noise. We assume that the dynamics of the detector lattice is well described by a three-dimensional 
distribution of acoustic phonons, with sound velocity v (the same for all three phonon modes). If the heat bath consists 
of a dielectric membrane, with the sound velocity w„i, then the phonon flux that penetrates through the separation 
surface, of area S*, from the detector into the membrane has the expression: 

(l^)'(^)'C(3)i, for 

Nk = { (26) 

^S^(^)'C(3)i, for«>t,„,, 

where t is the transmission coefficient, related to the acoustic impedances of the detector lattice, Z, and of the 
membrane, Z^, by the relation t — ^ZZ^J^Z + Z^^)^^ fl^l- The energy flux from the heat bath into the detector is 



47r2 vlh^ 



S^t-^{n ~ Tt) ^ ^KS{Ti - Tt), (27) 



and the power shot noise due to the quantisation of the phonon energy is 

(<5'QK,shot)^ = 8^Sk5(T| + Tf), (28) 

which, when Ti = T2, assumes the form {S^QK,shot)u^ = (C(5)/C(4))4fcBT2GK « AksT'^GK, where Gk = 
|9Qk/9Ti|ti=T2 is the Kapitza conductance. 
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c. Optical input power shot noise. Let us suppose that the detector admits a frequency band A/o from the 
electromagnetic spectrum within which the quantum efficiency is unity. If we denote by n{ujo) the density of phonons 
with the angular frequency ujo in a spatial region around the_ detector, then the energy flux and the number flux of 
photons on the detector are Qo{uJo) = ShLUoCn{uJo)Aujo/4:and No{uJo) = 5*071 (ijJo)AijJo/4, respectively, where ojo = 27r/o, 
as usual, c is the velocity of light, and we assume that AcJq <C lOq- The electromagnetic radiation interacts with the 
ES, therefore we shall take QojuJo) = Qoc('^o) and Qop(wo) = 0. Since n{uJo) — (l/7r^c^)a;o(exp [huJo/k^T) — 1)^^ (see 
for example Chap. 10 in Ref. Pq]), where T is the temperature of the background radiation, and if we assume that 
each photon produces a 5 peak of power into the ES, we obtain the following expressions for the input power and the 
shot noise into the ES of the TSE, respectively: 

n (::]-^ ^ ^° (90) 

^oo(UJo) - ^ ^2^2 exp (huJo/kBT) - 1 ^ ' 

and 

(<5^Qoc.hotK)>. = I f-) ' ,"°f"° , ■ (30) 

Z \TTC / exp (ftCJo/KBJ j — 1 

The radiation that influences the astrophysical measurements is the cosmic background radiation, which corresponds 
to T = Tb « 3 K. 

d. Shot noise in the junction. As anticipated in Section |, the shot noise fluctuation of the total power extracted 



through the NIS junction has no relevance to our problem. The calculations in Sections II A and II B showed in a 
rigorous manner that the quantity of interest is Qj, given by Eq. (^). The shot noise fluctuation of this quantity, 
which enters directly into the calculation of NEP, has the expression 

(<5'Qj,shot).; = ^ y deg{e)l{e~eV)^ 

X [/(e - eV, re)(l - 2/(e, T,)) + /(e, T,)] 

+ (e + eVflfie + eV, T,){1 - 2/(e, T,)) + /(e, r,)]| , (31) 

which is the term labeled by (^ in Section || The fluctuation of the particle flux, {S'^N ]^shot)uj, that enters the 
expression of T(uj), can be calculated in a similar way as {S'^ Q ,] .shot) ui , with the result 

2 f°° ( 

(^'^j,shot)c. = deg{e)Uf{e~eV,T,) + f{e + eV,T,)] 

x[l-2/(e,T3)] + 2/(e,T3)]|. (32) 

We have written down explicitly all the equations needed for the calculation of NEP. We are now left with the 
task of calculating the parameters of the refrigerating junctions for specific working conditions. Since the formulae 
involved are rather complicated P,p^ and the microbolometers work at low temperatures, we give below some useful 
analytical approximations, valid for this range of temperatures. 

III. APPROXIMATIVE ANALYTICAL EXPRESSIONS FOR THE CALCULATION OF THE 

JUNCTIONS PARAMETERS 

In order to calculate the efficiency of the microbolometer we have to know the parameters of the refrigerating 
junctions for given working conditions. The physical phenomena that take place in NIS microrefrigerators have been 



presented in a series of publications (see for example Refs. ]20|,|l|,|l0|, 2 11 19[|). Yet, the formulae involved are rather 
complicated and, in the general case, have to be calculated numerically. The metal that had been used extensively as 
a superconductor in the construction of the NIS refrigerating junctions is Al, which has an energy gap A « 200 /ieV 
at temperatures much lower than the critical one. Fortunately, the temperature range of interest for us (T < 0.3 K) 
is well below the critical temperature and since in this range A/k-^T ^ 1, we shall give here a general method to 
calculate analytical approximations for all the formulae presented in connection with the NIS junctions. 
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A. General formulae 



In the previous sections we saw that aU the formulae needed for the quantitative evaluation of the processes in NIS 
junctions [see Eqs. (|), (|), (i), (|), (0), d), (|l|), and (||)] have the general form: 

- /(e - E, T,))f{e, T,)] + r(e + E)' (33) 
X [/(e + E, re)(l - /(e, T,)) + - /(e + E, T,))f{e, T,)] 



e2i?T J A ^ %/e2 - A2 
(e - £;)'[/(£ - i^, Te) + s/(e, T,) - (1 + s)/(e - T,)fie, T,] 



X 

+r(e + EY [fie + E, T,) + sf{e, T,) - (1 + s)f{e + E, T,)f{e, T,)] 

The parameter / takes one of the values 0, 1, or 2, while r and s can be +1 or —1. The function 9{s) is the Heaviside 

step function. For example (+)F2^+' = (^^Qj ,hot)u., ^+^F^'^'' = (J^TVj^shot)^^, and = Qj. The integral in Eq. 

(||) splits up in an obvious way into six terms. We shall evaluate them one by one in the approximation 

A A 

A, = -—>A, = -—:»l. (35) 
fceJc kbJs 

Other notations that we shall use are: [3 = I/ZcbT, where any subscripts attached to T will be transferred to /3, 
5 = /\-E,a, = PJ, as = M: B, = /3e(A + E) :^ 1, and = /3s(A + E) » 1. 

The analytical approximations are possible due to the exponential dependence of the Fermi functions on energy, 
and are based on the Taylor expansion of the function 1/^/1 + x, around x = 0. If we also write (a + 6)' = a' + 
la'--'^b+ [l{l - l)/2ja'-262 (which is exact for I = 0, 1,2), we can calculate general analytical expressions for the six 
integrals in Eq. (|3J). In what follows we shall be interested just in the two highest order terms in A^, As, Be, and Bs 
in each of the integrals mentioned. For each /=0, 1, or 2 these highest order terms can be extracted rigorously from 
the general formulae given below, but in some cases the remaining lower order terms do not have the complete form 
(therefore we advise the reader to recalculate these lower orders, if they are necessary). 

Let us now start to calculate each of the integrals from Eq. (p4[). In these calculations, the values of the parameters 
r and s have no relevance, so they will be omitted. They should be introduced just in the end, when the expression 

of is to be reconstructed. The first term is 



- A2 exp [/3c (e - E)] + 1 

(fcBre)'+l 



/2a; 



{Ae [r(/ + l/2).g,+i/2(ae 



+Wil - l/2)aeg,_i/2(ac) + - 3/2)algi_yM 

+ ^ [T{1 + 3/2).g,+3/2(ae) + (Z + l/2)ae.g;+i/2(ac) 



ril-l/2)aigi_,f2{ae) 

where the functions gi{ct) are the Z*'' order polylogarithmic functions of argument — e~", gi{ct) 
[l/r(Z)] t'~^(e*+" + 1)~^ dt (see, for example, Refs. [gl|] and references therein for more details). 
In the second term: 

de. ^ 



A - A2 exp iPsc) + 1 

(fcBr.)'+i 



/2As 



{As [r{l + l/2)gi+,/2{As) 
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+ ^ [T{1 + 3/2).g,+3/2(^) + /r(; + l/2)a,g,+i/2(^) 



ra-l/2)a,Vi/2(As) 



since 5/(0;) = e " [1 + O (e ")], for any /, if a » 1, we shall use the approximation gi{As) 



-A, 



{A, [Til + 1/2) 



+W{1 - l/2)a, + t^r(; - 3/2)a2 



r(; + 3/2) + ir{i + i/2)as + ^^^^r^ - 1/2)0^ 



Within the same approximation, the third term, 



^3, 



de 



- A2 exp [/3e(e - E)] + I exp (/3,e) + 1 ' 



can be written in the form: 



h. - 2e- 



■E 

A;=0 



kl 



+W{1 + k~ l/2)ac.g;+fc-i/2(ac) + ^^^^r(; + k- ■i/2)algi+k-3/2iae) 
+ ^ [T{1 + k + 3/2).g,+fc+3/2(ac) + Zr(; + k + l/2)ae.g;+fc+i/2(ac) 



/(?-!) 



r(/ + /c-l/2)a,Vfe-i/2(ao) 



(36) 

If Tg 3> To we may just take the first term in the summation and obtain 

/3, « 2e~^^^h, . (37) 

This is a reasonable approximation, since for a typical case of Tc w 0.1 K and Tg ~ 0.4 K > T2 |jl^ w 6 and 
Ts/Tc ~ 4. In such a case /a, ^ /i, so there is no point in taking higher order terms in Eq. (|37|). If Tg = Tc, the infinite 
summations in Eq. (|36| ) do not converge. In such a case we can make use of the general inequality I^^ < 2e~''=^/ij, 
valid for any Tg, and eventually neglect this term, since in the cases of practical interest A^ > 20. 
The fourth term is 



de- 



(e + EY 



_ A2 exp [P,{e + E)] + 1 

(fcBTe)'+l 



/2X 



{A, [r{l + l/2)5,+i/2(i?e) 



+Wil - l/2)i?e<?,-l/2(i?e) + ^^^^r(Z - 3/2)B2g,_3/2(i?e) 
+ J [r(/ + 3/2).g,+3/2(i?e) + im + l/2)Be.g, + l/2(So) 



T{l~l/2)Blgi_,/^{B,) 



X A 



/2a: 

r(Z + 1/2) + IT{1 - 1/2)B, + t^r(Z - 3/2)^2 



r(; + 3/2) + zr(; + i/2)Be 



2 

/(?-i) 



r(; - 1/2)^2 
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The fifth term is 



dt 



eY 



Ve^ - A2 exp (/3se) + 1 



(fcBT.)'+i 

V2A 



{A, [ra + l/2).g,+i/2(A,) 



Aril-l/2)B,gi_,/M 
3 

'i 



ra-3/2)i?,V3/2(^) 



+ - [r(Z + 3/2).g,+3/2(A,) + Zr(; + l/2)B,gi+,/2(A,) 



'-Til-l/2)Btgi.u2{As) 



(fcBr.)'+^ 



{As [r{i + 1/2) 



+/r(z - i/2)Bs + 



1(1-1) 



Til - 3/2)B! 



+ - r{l + 3/2) + lT{l + l/2)Bs 
and the last term is 



1(1-1) 



Til-1/2)B! 



de 



(e + EY 



1 



Ve2 - A2 exp [f3eie + E)] 



(fcBTe)'+l 

V2a; 

r(/-3/2) 



r(r 



1 exp (/3s e 
1/2) 



+ 1- 



Til-1/2) 



(1 + Te/rs)'+i/2 (i + r,/Ts)'-V2 



+1 



2 (1 + Te/rs)'-V2 

r(; + i/2) ^ , ;g 



(1 + Te/Ts)' + l/2 



3 r r(/ + 3/2) 

4 [(l+Te/Ts)'+3/2 
ra - 1/2) ^2 



(i + re/rs)'-i/2 



One way to derive the last expression is to use the identity X]fc^o(~l)'^r7FnT'^'' ^ (^ 



^''WTT)''''' = (1 + c)"T(r), valid for any 
c G (—1, 1). For c = 1 the left hand side of the identity is not convergent, but wc can still define it as r(r)/2'', using 
the continuity property of the right hand side. 

Using the above analytical expressions, we can calculate easily the quantity (^)F;*-*-' for any I. As mentioned in the 
beginning of this section, in the concrete calculations we refer to the case in which the superconductor is Al, with 
A « 200 /zeV. Then for a temperature To = 0.1 K, Be « 46, so we shall neglect the terms I4, and /g, , since they 
are proportional to exp(— _Bo). Moreover, it should be noted that three of the other terms, I21, -^3, , and /s, , are all 
proportional to exp {—As). Among these three, /s, contains the highest orders for / = 1 and 2, while for Z = all of 
them are of the same order. Therefore care should be taken to calculate consistently the highest order of the terms 
proportional to exp (—As). 



B. Power flux through the NIS junctions 



In order to study the power flux through the NIS junctions, in the limit of low temperatures, let us take just the 
leading orders in Ae, As, i?c, and Bs from the formulae given for Ij^, j — 1, ... ,6 and I — 1. After dropping the 
irrelevant terms we arrive at the following expression ^ 



^ff3/2(ac) + acgi/2(ao) 




(38) 



A problem that occurs in all the cooling experiments by NIS junctions is the heating of the superconductor. Due to 
the poor heat transport properties of the superconductor, the temperature Ts increases too much and Q j becomes 
negative, or very close to zero, for any bias voltage. In such a case the cooler is practically ineffective. The heating 
can be controlled to some extent by placing normal metal traps on top of the superconductor. This has been discussed 
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in Refs. |^,|T9|. To avoid getting into too many details, here we shall consider Tg as an external parameter which is 
independent of the applied bias voltage (for any specific design more detailed calculations can be made). In case of 
an ideal quasiparticle trap, we set Tg — T2- For a constant To, Qj is just a function of a^. This function reaches 
its maximum value at = aopt.o ~ 0-66 [^|, which corresponds to what we shall call the optimum bias voltage, 
Vopt = A — A:Broaopt,o- From here we can draw an important conclusion: in the limit of low temperatures, the 
difference between the optimum bias voltage, Kjpt, and the energy gap A in the superconductor is independent of A 
and scales with the temperature of the ES in the normal metal. This means that 5opt = A — l^pt ^ as Tc ^ 0, and 
floptCTo) — > aopt.o as To 0, where aoptCTl;) is the exact value of at the maximum cooling power, at finite To. The 
ratio between Oopt and aopt,o &s a function of A^, is shown in Fig. |2[ We have to note that the formulae containing 
5fc(aopt.o) can be simplified by replacing this function neither by gfc(O) = |(1 — 2^~^)C,{k)\ nor by exp (— aopt,o) in the 
limit of low temperature. Replacing by aopt,o in Eq. (|38|), we obtain the optimum cooling power |jl|: 

A2 fk^T^V''^ A2 /27rfcBTg ^ , , 



^^opt/ ^opt,0 



. 9? 



.96 
. 94 - 
0. 92 
0.9- 



10 20 30 40 50e 

FIG. 2. The ratio between flopt, the value of at maximum cooling power and finite temperature Tc, and its limit value at 
zero temperature, fflopt.o, as a function of = A/ZcBle. Note that a range interesting for applications would be A^ > 20. 



IV. RESULTS 



Based on Eqs. (|2^) and (|19| ) and using the formulae given in Sections |II Q and|II[ we can make concrete calculations 
of NEP in the TSE of the microbolomctcr. In Ref. ||2^ the direct and indirect cooling methods were compared, for 
the case of zero bias and optical input power. The conclusion was that NEP was more than one order of magnitude 
lower in the case of indirect cooling than in the case of direct cooling, clearly recommending the first method for 
applications. When the optical input power is not zero the indirect cooling of the TSE might not be enough, due 
to the poor electron-phonon coupling at low temperatures [see Eq. (|2^)]. Therefore another set of junctions should 
be attached directly to the TSE to compensate for the optical power. To evaluate the situation, let us suppose 
that the resistances of these junctions are calibrated in such a way that Qj^opt = Qh + Qoe (we take Qop = 0), so 
To — Ti. We consider that Ti = r2 = 0.1 K and that the superconductor is in contact with an ideal quasiparticle trap 
(Tg = T2). The optical power incident on the detector is due to the cosmic background radiation and corresponds 
to a temperature Tb = 3 K. As seen form Eqs. (|9|) and (|3^), both Qod^o) and {5'^Qoc,shot{^o))u} depend on the 
detector dimensions, emissivity, bandwidth and frequency of the absorbed radiation. To eliminate the first three 
of these parameters, we shall use the formula ((5^0oo,shot(^o))a) = 2?icL;oOoc(wo) and treat Qoc('^oc) as a parameter 
obtainable from experiment. On the other hand, making use of the formulae given for the functions , P = 1, . . . , 6 
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and keeping only the highest orders in = A^, we obtain the following expression for the fluctuations: 

Now we can compare Eqs. ( ^9|) and (^). If we neglect the terms proportional to g-^/^BTc^ ^g^^ write in general 

(<5'Q,I,shot)^ ~ 3.47fcBTeQj,opt(^o) = 3.47fcBreQoo(Wo) . (41) 

Comparing the above equation with the equation for ((5^Qoc,shot(wo))(^, we obtain (i5^Qj,shot)c^ ~ 
1.74(fcBrc/?ia;)((5^(5oe,shot(wo))(^- If we write tOo = a^t^o,inax, where Wo,max is the angular frequency corresponding 
to the maximum of the energy density of the cosmic background radiation, ?ia;o,max = y^BTb {y ~ 2.82), then we can 
write 

(<5'Qj,shot)c. ~ O.G2^{6^Q,,,shotix))^ . (42) 



Since dQj/dV = at the optimum bias voltage, Eq. ( po[ ) simphfies and we obtain T(a;) = ((5^(5j,shot)w In experiments 
it is expected that Tb/Te > 30, so, for x of the order of 1 [Infra Red (IR) radiation] the shot noise in the junction, 
in the case of the combined cooling method, is very small as compared to the shot noise due to the background 
radiation. On the other hand it should be noted that the optical power of the IR radiation is absorbed by an antenna 
and is transformed into heat by the Joule effect in the TSE. In such a case, the noise contribution due to the input 



optical power may not be well described by Eq. (30) (here further investigation is needed). Therefore we shall write 
in general {S'^Qoc{^o))uj — '^(fi'huJoQoci^o), where is a parameter of value between and 1. 

Taking again the two limit situations from Ref. |22), and using Eq. ( ^ ) we say that if the Kapitza resistance is 
much smaller than the thermal resistance between the electrons and phonons, then 



NEP^ « (<52Qcp,shot)c. + (<5'QK,shot). 



0.624^) xykBThQoc{x) 

XlhJ 

j-2 A 1 ./,.„„„ Tc 



~ {S'Qep,shot)u + (^0 + 0-62;^ j xyfcBTbOoo(x), (43) 

where the last approximation holds if (TJj/Ti)^ is much smaller than the ratio between the electron-phonon and the 
Kapitza resistance [(T2/Ti)^(EnTi/EK5') <C 1], which is certainly the case for the indirect and combined cooling, 
under the assumptions made. In the other limiting situation, when the Kapitza resistance is much higher than the 
thermal resistance between the electrons and phonons, we can write 

NEP^ « 4(52Q^p^,hot)c. + (^0 + 0-62;^) xyksnQooix) . (44) 

Note that the coupling between the phonon and electron systems introduces a factor 4 in front of {S'^ Q cp, shot) lo in this 
situation. As in |2^] we assume that the volume of the TSE is fi = 1 /im'^ and if we take Ti = Tg = 0.1 K, we obtain 
(<5^Qop,shot)y^ ~ 7.43 X 10^^^ W/a/Hz. Let us suppose that we detect IR radiation and suppose further that in such 
a case 0^0. From Eqs. ( ^ and ( |4^ ) we can calculate two critical input optical powers, (5oo,i(a;) and Qoc,2{x), 
respectively, by the equation (5oc,i(2;) = Qoe,2{x) / A = ((5^Qep,shot)w/(0.62yfcB2l3). For our choice of temperatures we 
obtain Qoc,i(a;) « 2.29 x lO^^^ ^nd QocA^) ~ 9.15 x 10"" W. Therefore, for optical input power smaller than 
(500,1/2(2;) (depending whether we are in the first or in the second case), the noise contribution comes essentially 
from the electron-phonon coupling. For hig her input power, the noise due to the electron tunneling through the NIS 



junctions dominates. As seen from Eq. (42), in such a case the NEP has a weaker dependence on temperature, being 



1/2 

proportional to Tc , for constant optical input power. 

Throughout this paper we did not refer to the measurement of the temperature of the ES of the TSE. As mentioned 
in the beginning, there exist different methods to do that and we want to keep the arguments here as general as 
possible. Therefore, let us suppose that in an experiment the temperature is calculated from the measured value of a 
quantity Af (Tg) (for example M can be the voltage across the NIS junction biased at constant current, or the current 
in a voltage biased transition-edge sensor). If the inaccuracy in reading the quantity M is {5iM)^ (this should be 



12 



significant just in badly designed experiments) and the total mean square fluctuations of M (for example the shot 
noise added quadraticaly to the amplyfier noise) are ((5^Mf)^, then the uncertainity in the measurement of M is: 



{5Hl)]J^ = ((52 AW)^ + 2^Re({5M5T,)^) + {S^T,)^ + {5^M)^ , (45) 

where we assumed again that the changes in the chemical potential fi, due to temperature fluctuations, are negligible 
and Re({SM6Tc)uj) represent the real part of the correlations between the fluctuations of M and of Tq. To evaluate 
Eq. (^5|) we have to introduce in Eqs. (|l^) the noise contribution due to the bias power. In this way we obtain 
extra terms in the expressions for At{uj) given by Eqs. (|l^ ) and (|r^). Denoting the new function, influenced by 
the measurement process, by Amrriuj), and using in an obvious manner the two expressions for At(w), we write the 
equivalent of Eqs. (|l6|) and ( [l7| ) as 

A,„^(w) = 6T,iuj)ZNEP^ - At{lo) + ST.iiu)^ (46) 

and 

\ut{^) = At w + (5Qb,shot w - -TTT^-^ : 47) 

respectively, where Eh is the bias voltage of the thermometer and iVb.shot is the shot noise of the particle current due 
to the thermometer. The presence of the last term in Eq. ( |4^ ) depends on the temperature measurement setup. It 
should also be noted that when Qb is due to a flux of particles, so when is not zero, the 7r/2 angle between the 
phases of the noise spectra of these two quantities annihilates their correlations. Taking this into account, we can 
write in general 

(A^^(u;)) = {AUu)) + (52Q,,3hot). + ^ {^^^ {S'NKsi^oth 

In this situation, the total Noise Equivalent Power, NER, should be defined as the input power that would produce 
the same change in the quantity M as the one in Eq. (f45|). Simple calculations lead us to the general result 




-1 

Qj.^ \ZnepJ{SiM), 
where Znep^ was defined in Eq. (|46|). Using Eq. (|47|) we can write 



Re(^{5MfSTc)S^ ^ Re 



SMf 




For the case when the thermometer is a NIS junction, current or voltage biased, the NEP was calculated in Refs. 
@,[|j23| in the limit of zero Kapitza resistance, infinitely good reading accuracy of quantity M , and when the only 
refrigerating junction in the system is the thermometer junction. 



V. CONCLUSIONS 



In this paper we calculated the temperature fluctuations and the Noise Equivalent Power NEP in the thermal 
sensing element TSE of a microbolometer cooled by NIS junctions. From the general result we extracted simple 
expressions for two limiting cases. In the first case we considered that the Kapitza resistance is much smaller than the 
thermal resistance between the electrons and phonons, while in the second case we considered the opposite situation. 
In both cases the important role in the calculation of the NEP was played by the electron-phonon shot noise and the 
noise due to the cooling NIS junction, coupled directly to the TSE, and to the input optical power [see Eqs. ( ^ ) and 
(I)]- 
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We gave analytical expressions for the quantitative evaluation of the processes that take place in NTS junctions, 
which are valid in the range of low temperatures. The validity of these for mulae can be checked by the evaluation 
of the lower order terms in the quantities A^/c (see the definitions in Section III A| ). Using the analytical expressions 
we showed that in the limit of low temperature of the electron gas in the normal metal (Tc) the difference between 
the gap energy in the superconductor, A, and the bias voltage corresponding to the maximum cooling power, Vi^pt 
(supposing that the temperature of the superconductor does not vary with the bias voltage), is independent of A 
and scales with To, like A — Vopt = 0.66A:bTc- This is a very useful result since in such a case, in the quantitative 
evaluation of the processes that take place in NIS junctions, we can use neither the limit (A — Vopt)/kBTc 3> 1, nor 
the limit (A - Vopt) /ksT^ < 1, for -> 0. 

In the end we took into account in a general way the effect of the measurement process on the Noise Equivalent 
Power and we gave general expressions for its calculation. As expected, the measurement increases the uncertainty 
in the detected input optical power [see Eq. (Esh]. 
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